set more off

use "dataset_cities", clear

replace placepop = log(placepop + 1)

keep if sampled == 1

forvalues i = 50(50)200000 {
	local log`i' = log(`i' + 1)
}

#delimit;

hist 
	placepop
	,
		freq
		ytitle("Number of cities")
		ylab(, angle(horiz))
		xtitle("Population")
		plotregion(style(none))
		xlab(
			`log50' "50"
			`log100' "100"
			`log300' "300"
			`log1000' "1,000"
			`log10000' "10,000"
			`log50000' "50,000"
			`log200000' "200,000"
		)
		;
		
#delimit cr

gr export "figure A4 a.eps", replace
shell epstopdf "figure A4 a.eps"
